Submitted to ApJ 

Preprint typeset using I^T^X style cmulatcapj v. 8/13/10 



HIGH-REDSHIFT METALS. II. PROBING REIONIZATION GALAXIES WITH LOW-IONIZATION 

ABSORPTION LINES AT REDSHIFT SIX 1 ' 2 

George D. Becker 3 , Wallace L. W. Sargent 4 , Michael Rauch 5 , Alexander P. Calverley 3 

Submitted to ApJ 

ABSTRACT 

We present a survey for low- ionization metal absorption line systems towards 17 QSOs at redshifts 
z cm = 5.8 — 6.4. Nine of our objects were observed at high resolution with either Keck/HIRES 
or Magellan/MIKE, and the remainder at moderate resolution with Keck/ESI. The survey spans 
5.3 < z a bs < 6.4 and has a pathlength interval AX = 39.5, or Az = 8.0. In total we detect ten 
systems, five of which are new discoveries. The line-of-sight number density, £(X) — 0.25+q jjj (95% 
confidence), is consistent with the combined number density at z ~ 3 of DLAs and sub-DLAs, which 
comprise the main population of low-ionization systems at lower redshifts. This apparent lack of 
evolution may occur because low ionization systems are hosted by lower-mass halos at higher redshifts, 
or because the mean cross section of low-ionization gas at a given halo mass increases with redshift due 
to the higher densities and lower ionizing background. The roughly constant number density notably 
contrasts with the sharp decline at z > 5.3 in the number density of highly- ionized systems traced by 
C IV. The low-ionization systems at z ~ 6 span a similar range of velocity widths as lower-redshift 
sub-DLAs but have significantly weaker lines at a given width. This implies that the mass-metallicity 
relation of the host galaxies evolves towards lower metallicities at higher redshifts. These systems 
lack strong Si IV and C IV, which are common among lower-redshift DLAs and sub-DLAs. This is 
consistent, however, with a similar decrease in the metallicity of the low- and high-ionization phases, 
and does not necessarily indicate a lack of nearby, highly-ionized gas. The high number density of 
low-ionization systems at z ~ 6 suggests that we may be detecting galaxies below the current limits 
of i-dropout and Lya emission galaxy surveys. These systems may therefore be the first direct probes 
of the 'typical' galaxies responsible for hydrogen reionization. 

Subject headings: cosmology: observations — cosmology: early universe — intergalactic medium — 
quasars: absorption lines 



1. INTRODUCTION 

The universe at redshift six remains one of the most 
challenging observational frontiers. The high opacity of 
the Lya forest inhibits detailed measurements of the in- 
tergalactic medium (IGM), which hosts the vast major- 
ity of baryons in the early universe. At the same time, 
extreme luminosity distances make it difficult to assem- 
ble large samples of galaxies, or to study individual ob- 
jects in detail. Despite these challenges, however, sig- 
nificant observational progress is being made at z ~ 6. 
Statistics of the scarce transmitted flux in the Lya for- 
est are being combined with insights from quasar prox- 
imity zones_to learn about the very high-redshift IGM 
(e.g.JFan et al.|l2006tlLidz et alJl200alBecker et al.ll2007b 
Gallerani et alJ 120081: iMaselli et all l2009t iBolton et all 
20101 : iCalverlev et al.ll2010HCarilli et al.ll201fJ : iMesingerl 
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fornia Institute of Technology and the University of California; 
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Chile. 

3 Kavli Institute for Cosmology and Institute of As- 
tronomy, Madingley Road, Cambridge, CB3 0HA, UK; 
gdb@ast.cam.ac.uk, acalver@ast.cam.ac.uk 

4 Palomar Observatory, California Institute of Technology, 
Pasadena, CA 91125, USA; wws@astro.caltech.edu 
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2010). Deep galaxy surveys from the ground and from 
space are also beginning to yield considerable informa- 
tion o n the properties of g a laxies at z ~ 6 and beyond 
(e.g.. IBouwens et all 120071 120101: IStanway et ail 120071: 
Martin et al.H2008tlOuchi et al.ir2 008. 20 09HBunker et al.1 
20101: IHu et all I2010t ILabbe et all 2010t iMc Lure et al.1 
2010: IQesch etail 2010: St ark et al.ll2010f) . A more com- 
plete understanding of the universe at such early times, 
however, will require insight from a wide variety of ob- 
servations. 

Metal absorption lines offer a unique probe into the 
high-redshift universe. Transitions with rest wave- 
lengths longer than H I Lya will remain visible in QSO 
and GRB spectra even when the forest becomes sat- 
urated. While absorption lines studies at z ~ 6 are 
still in their early days, metal lines are expected to 
provide insight into past and ongoing star formation, 
the action of galactic winds, as well as conditions in 
the interstellar medi a of galaxies, mst as they do at 
lower redshifts (e.g . . lSongailall200ll; Simcoe et al.l l2002t 
Schave et all [20031 lAdelberger et al 1 120051: iWolfe et al.l 
20051: IScannapieco et all 120061 : iFox et al.ll2007D . At high 
redshifts, the role of metal absorption lines in study- 
ing galaxy evolution becomes increasingly significant as 
galaxies become more difficult to study in emission. 

Metal lines play another unique role at z > 6 as 
potential probes of hydrogen reionization. In a sce- 
nario where galaxies provide most of the ionizing pho- 
tons, dense regions of the IGM may become chemically 



2 



Becker et al. 



enriched early on, yet remain significantly neutral un- 
til the end of reionization due to the short recombi- 
nation times at high densities. High-density peaks are 
commonly predicted to host the last pockets of neu- 
tral gas following the growth and o v erlap of large H II 
regions jM iralda- Escude et al.1 l2000t iCiardi et all 120031 : 
iFurlanetto fc Qhll2005t iGnedin fc Fanll2006D . If these re- 
gions are metal-enriched, then forests of low-ionization 
absorption lines such as O I, Si II, and C II m ay be ob- 
servable during this last phase of reionization (|Ohll2002l ; 
IFurlanetto fc Loebl[200l . 

We have conducted a two-part survey for metals out 
to z ~ 6. In the hrst paper of this series (|Becker et al.l 
120091 hereafter Paper I) we used high-resolution, near- 
infrared spectra to search for C IV over z = 5.3 — 6.0. 
C IV is a widely-used t racer of highly -ionized metals in 
the IGM at z < 5 (e.g.. lSonga"ilall2001| ). and initial stud- 
ies had ind icated that C IV may remain abundan t out 
to z ~ 6 (|Simcoel [200l iRvan- Weber et al.l [20061) . In 
contrast, the longer path length and greater sensitivity 
of our data allowed us to demonstrate that the number 
density of C IV systems declines by at least a factor of 
four fro m z ~ 4.5 to z > 5 .3. T his downturn was also 
seen by IRvan- Weber et al.l |2009) using a larger sample 
of low-resolution spectra. Since C IV is expected to be 
an effective tracer of metals in the IGM at z ~ 5 — 6 
( Oppenheimcr fc Pavel [200l lOppenheimer et al J 12009), 
the decline suggests that there are fewer metals in the 
IGM towards higher redshifts, consistent with gradual 
enrichment of the IGM by star-forming galaxies. The 
apparent rapid evolution of C IV near z ~ 5 may also 
point to changes in the ionization state of the metal- 
enriched gas, perhaps due to a hardening of the UV back- 
grou nd associated with the b eginning of He II reioniza- 
tion (jMadau fc Haardtll2009T ). 

In this paper we present a complementary search for 
low-ionization metal absorbers over 5.3 < z < 6.4. Thi s 
is an extension of an earlier work (IBecker et al.l [2006). 
in which we used high-resolutions spectra of nine QSOs 
spanning 4.9 < z cm < 6.4, including five at z cm > 5.8, to 
search for 'Ol absorbers,' so named due to the presence 
of O I A1302. Although we identified a possible excess 
of these systems at z > 6, the significance of that re- 
sult remained unclear due to the relatively short survey 
pathlength, and the fact that the z ~ 6 systems occurred 
almost entirely along a single sightline. Here we increase 
our sample of z em ~ 6 QSOs to 17. The expanded survey 
offers us the opportunity to gather a more representative 
sample of low-ionization systems, providing insight into 
early galaxies and possibly the tail end of hydrogen reion- 
ization. 

The remainder of the paper is organized as follows. 
The data are described in Section [2j We present our de- 
tections and the measurements for individual systems in 
Section [3J In Section 2] we analyze a number of ensemble 
properties, including the number density, ionic mass den- 
sities, relative abundance patterns, lack of strong high- 
ionization lines, and velocity width distribution. We also 
present evidence that the mass-metallicity relation of the 
host galaxies of low-ionization absorbers evolves out to 
z ~ 6. Possible causes for the apparent lack of evolu- 
tion in the number density of low-ionization absorbers 
are discussed in Section where we also consider the 
likely hosts of these systems. Finally, we summarize our 



conclusions in Section [51 Throughout th is paper we as- 
sume (Q m ,n A ,h) = (0.274,0.726,0.705) ([Komatsu et al.l 
[20091) . 

2. THE DATA 

We have obtained high- or moderate-resolution spec- 
tra of 17 QSOs with emission redshifts z cm = 5.8 — 6.4. 
This sample contains more than three times the num- 
ber of z > 5 . 8 QS O sight lines than were included in 
IBecker et al.l (J2006) . Nine of our objects (generally the 
brightest) were observed at high resolution with either 
the HIRES spectrograph dVogt et al.[|1994l) on Keck or 
the MIKE spectrograph (|Bernstein et al.ll2002i ) on Mag- 
ellan. This sample includes the five z > 5.8 QSOs from 
IBecker e~a l. (2006), although additional data have since 
been acquired for some of these. We complement our 
high resolution sample with a set of eig ht QSO spec- 
tra ta ken at moderate resolution with ESI ([Sheinis et al.l 
2002) on Keck. The data are summarized in Table [1] 

All data were reduced using custom sets of IDL rou- 
tines that inc luded optimal sky subtraction algorithms 
([Kelson! [20031 ). In order to efficiently reject cosmic rays 
and bad pixels, a sin gle 1-D spectrum was optimally ex- 
tracted (|Hornelll986l) from all 2-D sky-subtracted frames 
for a given object simultaneously. Prior to extraction, 
relative flux calibration across the orders was performed 
using response functions generated from standard stars. 
Continua redward of the Lya forest were fit by hand us- 
ing a cubic spline. 

The HIRES spectra were taken using the upgraded de- 
tector. Two grating settings with the red cross disperser 
provided continuous spectral coverage up to nearly 1 /im. 
The spectra were taken using a 0'/86 slit, which provides 
a resolution FWHM of 6.7 km s -1 . The final ID spectra 
were binned using 2.1km s _1 pixels. Telluric absorption 
corrections based on a standard star were also applied 
to the individual 2-D frames prior to final extraction. 
For objects observed at multiple times during the year, 
changes in heliocentric velocity allowed telluric features 
to b e efficiently exclude d from the final spectra. Data 
from IBecker et all ([20061 ) were re-reduced to take advan- 
tage of the improvements in the reduction pipeline. Only 
data taken with the upgraded detector were included. 

MIKE spectra were taken using a single setting that 
gave continuous coverage up to 9400 A. We used a l'.'OO 
slit, which gives 13.6 km s _1 resolution in the red. The 
ID spectra were binned using 5.0 km s _1 pixels. Tel- 
luric corrections were performed after extracting the 
ID spectra using a high-resolution, high signal-to-noise 
transmission spectr um derived from solar observations 
([Hinkle et al.1 12003). which was shifted, smoothed, and 
scaled in optical depth to match the QSO spectra. 

ESI provides complete wavelength coverage up to 1 fj,m 
in a single setting. We used both the 0'.'75 and l'.'OO slits, 
depending on the seeing. Our combined resolution was 
50 km s _1 , and was roughly consistent across all objects. 
The ID spectra were binned using 15.0 km s" 1 pixels. 
Telluric corrections were performed after extraction us - 
ing the transmission spectrum from lHinkle et all {2003). 

Finally, we include Keck NIRSPEC echelle spec- 
tra (FWHM = 23 km s _1 ) covering high-ionization 
lines (Si IV and C IV) for SDSS J0818+1722 and 
SDSS J1148+5251. These data are described in Paper I. 
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Table 1 

Data 



QSO 


7 a 


Instrument 


^cxp 

(hrs) 


S/N tes h 




Hif 


;h Resolution 








nncc iac\aa mot; 

bJJbb J1U44— U125 


5.782 


MIKE 


7.5 


10.7 - 


19.7 


bJJbb JUoob+UUo4 


5.81 


HIRES 


12.5 


15.0 - 


23.4 


bDbb J0002+2550 


5.82 


HIRES 


14.2 


13.8 - 


28.0 


bDbb J0818+1722 


6.00 


HIRES 


8.3 


11.8 - 


18.5 


bDbb J 1306+0356 


6.016 


MIKE 


6.7 


7.3 - 


15.9 


bDbb J1048+4o37 


6.23 


HIRES 


5.0 


7.2 - 


12.9 


SDSS J1623+3112 


6.247 


HIRES 


12.5 


6.0 - 


20.6 


SDSS J1030+0524 


6.308 


HIRES 


10.0 


6.7 - 


16.9 


SDSS J1148+5251 


6.419 


HIRES 


14.2 


10.5 - 


21.6 




Moderate Resolution 








SDSS J0005-0006 


5.85 


ESI 


1.2 


27.6 - 


50.0 


SDSS J0203+0012 


5.85 


ESI 


1.8 


14.0 - 


22.8 


ULAS J0148+0600 


6.00 


ESI 


3.1 


50.9 - 


85.9 


SDSS J1630+4012 


6.05 


ESI 


0.6 


14.7 - 


27.7 


SDSS J0353+0104 


6.05 


ESI 


3.7 


27.8 - 


55.4 


SDSS J2054-0005 


6.06 


ESI 


3.3 


33.3 - 


50.8 


SDSS J2315-0023 


6.12 


ESI 


6.7 


27.5 - 


70.6 


CFHQS J0050+3445 


6.25 


ESI 


1.7 


21.7 - 


55.4 



a QSO redshifts q uoted to three decimal places are Mg II or CO redshifts 
from Carilli et al. (2010). Others are discovery redshifts based on Lya+N V 
emission, or are measured from the apparent start of the Lya forest. 
b Interquartile values of signal-to-noise in the continuum per resolution element 
within the wavelength interval used to search for O I. The resolution elements 
for HIRES, MIKE, and ESI are 6.7, 13.6, and 50 km s _1 , respectively. 

3. DETECTED SYSTEMS 

At z < 5, low- ionization absorption systems are typi- 
cally identified by their strong H I absorption, particu- 
larly the extende d damping wings of damped Lya sys- 
tems (DLAs; e.g- JWolfe et al.ll2005D . At higher redshifts, 
however, absorption features in the Lya forest become 
sufficiently blended that identifying individual absorbers 
becomes difficult. To identify low-ionization systems, we 
therefore searched redward of the forest for groups of 
metal lines including Si II A1260, Si II A1304, O I A1302, 
and C II A1334. A detection required that two or more 
lines occur at the same redshift. The lines were required 
to have similar velocity profiles, but the relative strengths 
of the different elements were not constrained. For each 
sightline we searched between the QSO redshift and the 
redshift at which O I enters the Lya forest. The mini- 
mum and maximum redshifts covered are 5.33 and 6.42. 

In total we detected ten sys tems, five of which were 
reported in iBecker et al.1 (|2006[ ) . Each contains Si II and 
C II, and all but one contains detectable O I. Seven of the 
systems occur along sight lines for which we have high- 
resolution data. The remainder fall towards QSOs we 
observed with ESI. Five of the seve n systems in the hig h 
resolution sample were reported in IBecker et al.l ((2006), 
while the two additional high-resolution systems occur 
along a single sightline, S PSS J0818+1722 Thes e two 
systems were also noted by iD'Odorico et al.l (|201lD . The 
three systems in the ESI data occur along separate sight 
lines. 

The systems are shown in Figures [TlfTOl In addition 
to the low-ionization lines (Si II, O I, and C II), we plot 
the wavelength regions covering Si IV and C IV in cases 
where we have data. Si IV may fall either in the optical or 
near-infrared, depending on the redshift of the system. 
For C IV, our coverage comes entirely fr om our NIR- 
SPEC data (Paper I). The systems from IBecker et al.l 
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Figure 1. Stacked velocity plot for the absorption systems at 
z = 5.7911 towards SDSS J0818+1722. Data covering the low- 
ionization lines and Si IV are from HIRES, while the data for C IV 
are from NIRSPEC. Histograms show the normalized flux, and 
the formal error array is plotted along the bottom of each panel. 
Shaded regions indicate detected transitions. 
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Figure 2. Stacked velocity plot for the absorption systems at 
z = 5.8765 towards SDSS J0818+1722. Data covering the low- 
ionization lines and Si IV are from HIRES, while the data for C IV 
arc from NIRSPEC. Lines are as in Figure Q] Si II A1260 is partially 
affected by a skyline residual. 




-100 

A v (km s -1 ) 



100 



Figure 3. Stacked velocity plot for th e absorption syste ms at 
z = 5.8415 towards SDSS J1623+3112 l|Becker et al.ll200r3 V All 
data are from HIRES, although the the data covering Si IV are 
binned to 8.4 km s — 1 pixels for display. Lines are as in Figure [T] 
C II in the bluest component is blended with a strong Mg I line at 
lower redshift. O I and Si II A1304 are both partially blended with 
a C IV system at z = 4.754, which is visible over —50 kms -1 < 
Av < kms -1 . 
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Figure 4. Stacked velocity plot for t he absorption syst ems at 
z = 6.0115 towards SDSS J1148+5251 (Becker et al. 2006). Data 
covering the low-ionization lines are from HIRES, while the data 
for Si IV and C IV are from NIRSPEC. Lines are as in Figure [T] 
The absorption covering the blue half of the C IV A1548 region is 
unrelated Fc II at lower redshift. 



( 2006) are reproduced primarily to demonstrate the ap- 
parent lack of high-ionization absorption lines, a point 
to which we will return later. In the next sections we 
describe our basic measurements and completeness esti- 
mates. 

3.1. Equivalent Widths 

We measured the rest-frame equivalent widths, Wo, for 
each available transition. The results are given in Ta- 
ble [2] For detected low-ionization lines, the wavelength 
intervals over which the equivalent widths were measured 
are shown as shaded regions in Figures [D llOl For non- 
detected low-ionization lines, upper limits are measured 
over the same intervals. We define our upper limits to be 
the sum of any net absorption, as in the case of blends 
with unrelated lines, and the uncertainty in the normal- 
ized flux integrated over the chosen window. 

High-ionization lines are often kinematically distinct 
from low-io nization lines in l ower- redshift DLAs and sub- 
DLAs fe.g.. lFox et~aLll2007t ). We therefore place conser- 
vative upper limits on the equivalent widths of Si IV 
and C IV by integrating over a ±100 km s _1 window 
around the nominal system redshift. We measure limits 
for both lines in the doublet separately, but take ad- 
vantage of the fact that Wo for the stronger transition 
should be no more than twice that of the weaker tran- 
sition. For Si IV, the smaller of Wo(Si IV A1394) and 
2 W (Si IV A1403) is given in Table [2 For C IV, we give 
the smaller of W (C IVA1548) and 2W (C IVA1551). 
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Table 2 

Equivalent Widths 





QSO 




Si II A126U 


Si II A1304 


(J I A13U2 


/" ' tt \ 1 Q 1 A 

O II A1334 


bi IV A1394 


Vj IV A1548 


SDSS 


J1630+4012 


5.8865 


0.235 


< 0.076 


< 0.031 


0.240 


< 0.362 




SDSS 


J2054-0005 


5.9776 


0.198 


0.072 


0.124 


0.221 


< 0.096 




SDSS 


J2315-0023 


5.7529 




< 0.195 


0.238 


0.240 


< 0.069 




SDSS 


J0818+1722 


5.7911 


0.267 


0.056 


0.182 


0.230 


< 0.105 


< 0.072 


SDSS 


J0818+1722 


5.8765 


0.064 


< 0.008 


0.058 


0.056 


< 0.044 


< 0.038 


SDSS 


J1623+3112 


5.8415 




0.135 


0.391 


0.278 


< 0.060 




SDSS 


J1148+5251 


6.0115 




0.038 


0.162 


0.154 


< 0.063 


< 0.160 


SDSS 


J1148+5251 


6.1312 




0.016 


0.079 


0.067 


< 0.015 


< 0.040 


SDSS 


J1148+5251 


6.1988 


0.020 


< 0.028 


0.020 


0.015 


< 0.072 


< 0.H3 


SDSS 


J1148+5251 


6.2575 


0.045 


< 0.008 


0.036 


0.047 


< 0.030 





Note. — Rcst-framc equivalent widths are given in Angstroms. Upper limits are lcr. 
a Minor differences from redshifts given in Becker et al. (2006) are due to an improvement in the HIRES wavelength 
solution. 

b This gives the smaller of the upper limit on Si IV A1394 and twice the upper limit on Si IV A1403. 
c This gives the smaller of the upper limit on C IV A1548 and twice the upper limit on C IV A1551. 
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Figure 5. Stacked velocity plot for t he absorption syst ems at 
z = 6.1312 towards SDSS J1148+5251 IjBecker et all 20061. Data 
covering the low-ionization lines are from HIRES, while the data 
for Si IV and C IV are from NIRSPEC. Lines are as in Figure [T] 



1.0 
0.5 
0.0 
1.0 
0.5 
0.0 
1.0 
0.5 
0.0 

1.0 
0.5 
0.0 
1.0 
0.5 
0.0 

1.0 
0.5 
0.0 





'•Or 

0.5 
0.0 



Si IV 1393.76 



C U M 334.53 




Sill Jl 260.42 



-100 100 

A v (km s -1 ) 

Figure 6. Stacked velocity plot for t he absorption syst ems at 
z = 6.1988 towards SDSS J1148+5251 (Bicker et al. 2006). Data 
covering the low-ionization lines are from HIRES, while the data 
for Si IV and C IV are from NIRSPEC. Lines are as in Figure [l] 



3.2. Column Densities 

Column densities of each ion were measured for sys- 
tems in the high resolution sample, where the velocity 
structure is most likely to be resolved. For six out of 
seven of these, column densities w ere measured using 
the a pparent optical depth method (Sa vage fc Sembachl 
1991). Here the optical depths of individual pixels, 7V, 
are measured, and the total column density in cm~ 2 is 
given by 



N = 



3.768 x 10 14 



/A 



(1) 



where / is the oscillator strength, Ao is the rest wave- 
length in Angstroms, and Svi is the velocity width of a 
pixel in km s _1 . This approach has the advantage that 
is does not require the absorption to be decomposed into 
individual components. We integrated the optical depths 
over the shaded regions in Figures [T][7] Prior to measur- 
ing the low-ionization lines, we smoothed our spectra by 
convolving the flux with a Gaussian kernel whose FWHM 
is equal to one half of the instrumental resolution. This 
takes advantage of the fact that the spectra are known to 
have a finite resolution, which is roughly three pixels for 
our choice of binning. Mild smoothing prevents individ- 
ual noisy pixels with flux near or below zero from dom- 
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Table 3 

Column Densities for High-Resolution Systems 



QSO 




log N Si „ 


log N Q 1 


log N c „ 


log N Si IV 


log iVc iv 


SDSS J0818+1722 


5.7911 


13.41 ± 0.03 


14.54 ±0.03 


14.19 ±0.03 


< 12.9 


< 13.4 


SDSS J0818+1722 


5.8765 


12.78 ±0.05 


14.04 ±0.05 


13.62 ±0.07 


< 12.6 


< 13.0 


SDSS J1623+3112 


5.8415 


14.09 ±0.02 


> 15.0 


> 14.3 


< 13.0 




SDSS J1148+5251 


6.0115 


13.51 ±0.03 


14.65 ±0.02 


14.14 ±0.06 


< 12.7 


< 13.8 


SDSS J1148+5251 


6.1312 


13.29 ±0.09 


14.79 ±0.34 


13.88 ±0.22 


< 12.3 


< 13.2 


SDSS J1148+5251 


6.1988 


12.12 ± 0.05 


13.49 ±0.13 


12.90 ±0.11 


< 12.9 


< 13.6 


SDSS J1148+5251 


6.2575 


12.67 ± 0.03 


13.97 ±0.15 


13.56 ±0.14 


< 12.5 





Note. — Column densities were computed using the apparent optical depth method in all cases except the 
z = 6.1312 system towards SDSS J1148+5251. For this system, column densities were measured using Voigt profile 
fits, for which a Doppler parameter b = 5.7 ± 1.0 km s — 1 was returned. 




-100 100 

A v (km s"') 

Figure 7. Stacked velocity plot for t he absorption syst ems at 
z = 6.2575 towards SDSS J1148±5251 jBecker et aL 2006). Data 
covering the low-ionization lines are from HIRES, while the data 
for Si IV are from NIRSPEC. Lines are as in Figure Q] C II is 
partially affected by skyline residuals. 



mating the summed optical depths, while only slightly 
degrading the spectral resolution (here by ^12%). In 
cases where the spectra had high S/N, or for transitions 
which contained no pixels near zero, smoothing had a 
negligible effect on the measured column density. We 
note, however, that smoothing should be used cautiously 
where the flux may genuinely go to zero. For example, 
in the z = 6.1312 absorber towards SDSS J1148+5251, 
O I A1302 is marginally saturated. We therefore fit Voigt 
profiles rather than apply the apparent optical depth 
method in this case. This system appears to have a single 
narrow component, and so by requiring that each transi- 
tion have the same redshift and Doppler parameter (i.e., 
assuming turbulent broadening), we can obtain a rough 
column density estimate for O I even though the central 
two pixels show nearly zero flux. 

Upper limits on the column densities were computed 
for non-detected species. For low-ionization lines, the 
optical depths were integrated over the same velocity in- 
terval where other low-ionization species were detected. 




-200 -100 100 
A v (km s"') 



200 



Figure 8. Stacked velocity plot for the absorption systems at 
z = 5.8865 towards SDSS J1630±4012. The data are from ESI. 
Lines a re as in Figure [l] Note the change in velocity scale from 
Figures [ll7l 

Limits on C IV and Si IV were set by integrating the 
apparent optical depths over an interval -100 km s _1 to 
+ 100 km s _1 from the optical depth-weighted mean red- 
shift of the low-ionization lines. For the high-ionization 
doublets, the column density in the zth velocity bin was 
taken to be the inverse variance-weighted mean of the 
column densities measured from both transitions in cases 
where the measurements were consistent within the er- 
rors. In cases where they were not consistent, as would 
be expected if one transition was blended with an unre- 
lated line, the smaller of the two values was taken. 

We used our optical depth measurements of Si II to 
infer another quantity used to classify low-ionization 
systems, the equ i valent width of Si II A 1526, W1526 
( Prochask a et al.l [2008) . Since we do not cover 
Si II A1526 with our current data, the lines profile for 
this transition was reconstructed from the optical depths 
of Si II A1260 and/or A1304. In all cases where this was 
done, at least one of these transitions was covered and 
not saturated. Since the optical depth should scale lin- 
early with the oscillator strength, we therefore expect 
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Figure 9. Stacked velocity plot for the absorption systems at 
z = 5.9776 towards SDSS J2054-0005. The data are from ESI. 
Lines a re as in Figure [T] Note the change in velocity scale from 
Figures 11171 
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-200 -100 100 200 
A v (km s -1 ) 

Figure 10. Stacked velocity plot for the absorption systems at 
z = 5.7529 towards SDSS J2315-0023. The data are from ESI. 
Lines a re as in Figure [T] Note the change in velocity scale from 
Figures IH7I Si II A1304 is blended with an unrelated absorption 
line. 

our estimates of W1526 to be reasonably accurate in all 
cases. The results are given in Table 01 

3.3. Velocity Widths 

Finally, we measured the total velocity width of sys- 
tems in our high resolution sample. We use the veloc- 
ity interval covering 90% of the optical depth in low- 
ionization species, Augo, which has been used to charac- 
terize the kinematics of DLAs and sub-DLAs at lower 
redshifts (jProchaska fc Wolfd I1997D . The widths are 
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13.5 14.0 
log N(C II) 
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Figure 11. Completeness estimates as a function of C II column 
density. The solid line is the pathlength-weighted mean for the 
full the sample, while the dashed and dotted lines arc for the high- 
resolution and ESI data, respectively. See text for details. 

measured from the shaded pixels in Figures [TJT] Mild 
smoothing was again applied. The results are summa- 
rized in Table 3) The transition used to measure Atig 
were chosen to have significant optical depth without be- 
ing obviously saturated. O I was used only in cases where 
the velocity widths of O I, C II, and Si II are apparently 
equal. In other cases, C II and Si II appear to have ab- 
sorption at high velocities not present in O I, presumably 
due to a change in the ionization state of the gas. 

3.4. Completeness 

We estimated our completeness by attempting to re- 
cover artificial absorption systems inserted into the data. 
The synthetic lines were modeled as Voigt profiles with 
a Doppler parameter of 10km s _1 (At>go = 23km s _1 ). 
This choice focuses on the detectability of the narrower, 
weaker lines in our sample. Artificial systems containing 
O I, C II, and Si II were placed at random into the data, 
assuming column density ratios equal to the mean val- 
ues measured in the data (Section FO]) . We then used an 
automatic routine to determine whether these systems 
would be detected. The detection criteria were tuned to 
mimic the results of a by-eye identification. A detection 
first required that a peak flux decrement greater than 
lcr be present at the same redshift for at least two of 
Si II A1260, O I A1302, and C II A1334 after the flux 
was smoothed by convolving with a Gaussian kernel of 
FWHM equal to the instrumental resolution. For can- 
didate pairs, the unsmoothed data were then fit using 
Gaussian profiles. A detection required that the ampli- 
tude and FWHM for both lines were measured at greater 
than 3cr confidence, and that the measured redshifts and 
velocity widths agreed with high confidence. The ratio 
of the maximum optical depths in the two lines were also 
required to be within a factor of 0.4 to 2.0 of the ratio 
expected from the assumed relative abundances. 

Each line of sight was tested with 10 4 artificial sys- 
tems over a range of column densities and covering the 
entire redshift interval used to search for low-ionization 
systems. The pathlength-weighted mean results for the 
full sample are shown in Figure 1111 along with separate 



Becker et al. 



Table 4 



Line Widths for High-Resolution 


Systems 




QSO 


-^abs 


Au 9 o a 


Transition'' 


Wl526° 






(km s — 1 ) 




(A) 


SDSS J0818+1722 


5.7911 


134 


C II A1334 


0.066 


SDSS J0818+1722 


5.8765 


25 


O I A1302 


0.016 


SDSS J1623+3112 


5.8415 


223 


Si II A1304 


0.249 


SDSS J1148+5251 


6.0115 


41 


O I A1302 


0.074 


SDSS J1148+5251 


6.1312 


17 d 


C II A1334 


0.028 


SDSS J1148+5251 


6.1988 


17 


Si II A1260 


0.004 


SDSS J1148+5251 


6.2575 


13 


Si II A1260 


0.012 



a Velocity width spanning 90% of the optical depth in the low-ionization lines 
b Transition used to measure Avqq 

c Rest equivalent width of Si II A1526 estimated from optical depth measure- 
ments of available Si II lines 



d Measured from Voigt profile fit 



results for the high-resolution and ESI data. We estimate 
that we are 50% complete at logiV(C II) ~ 13.3 in the 
high-resolution data and logiV(C II) ~ 13.7 in the ESI 
data. Note that the completeness estimates do not go to 
100% for large column densities due to the fact that parts 
of the spectra are unusable due to strong atmospheric ab- 
sorption or the presence of other absorption lines. This 
is particularly an issue at redshifts where Si II A1260 
is in the Lya forest, and the system must be identified 
from O I A1302 and C II A1334 alone. Realistically, how- 
ever, strong systems tend to have multiple components, 
making them easier to detect, and often have significant 
Si II A1304 absorption. 

4. SAMPLE PROPERTIES 

4.1. Number Density 

We now address the number density of low-ionization 
absorption systems at z ~ 6. We compute the number 
density of absorbers per unit absorption pathlength in- 
terval, t(X) = dAf/dX. Here, X is defined as 



dz . 



(2) 



where H(z) is the Hubb le constant at a given redshift 
(|Bahcall fc Peebles! Il969h . A non-evolving population of 
absorbers will have a constant number density per unit 
X. Our total pathlength is AX = 39.5 (Az = 8.0, 
or Al = 3.6 comoving Gpc). The high-resolution data 
cover AX = 21.0, while the ESI data cover AX = 18.5. 
Uncorrected for completeness, the line of sight number 
density in the entire sample is £(X) = 0.251^3 (95% 
confidence), where the error bars include Poisson uncer- 
tainties only. 

Corrections to £(X) for completeness are small com- 
pared to the errors, except when considering our 
weakest system. The z = 6.1988 system towards 
SDSS J1148+5251, which has logiV(C II) = 12.9, would 
have been detected along only 9% of our pathlength. Few 
systems that are similar are known at lower redshifts. 
The only comparable system of w hich we are aware i s 
the sub-DLA at z = 3.7 found bv IPeroux et all (|2007f ). 
which has logiV(C II) = 13.1 yet still contains O I. It is 
intriguing to consider whether more weak systems may 
be found at z ~ 6 in future, more sensitive surveys. 

Seven out of our ten systems are detected in the high- 
resolution data, which suggests that the two lines of 



sight with multiple systems may be unusual. At a given 
signal-to-noise ratio, however, HIRES and MIKE are 
more sensitive to narrow absorption systems than ESI, 
which has considerably lower resolution. In order to de- 
termine whether the two data sets are consistent, we 
first take the ratio of the completeness estimates in Fig- 
ure [Tl] at the column densities of the seven systems in our 
high resolution sample. Significantly, the weakest sys- 
tem (logJV(C II) = 12.90) would be nearly undetectable 
in the ESI data, while the completeness of the ESI data 
for the two systems with log N (C II) ~ 13.6 would be 
roughly half that of the high-resolution data. If we also 
take into account the slightly shorter pathlength of the 
ESI data, and the fact that the ESI data are somewhat 
less sensitive to even the stronger systems, then we would 
expect to detect ^4 systems in the ESI data if the high- 
resolution systems represent a fair sample. The fact that 
we detected three systems with ESI suggests that the two 
data sets are broadly consistent. 

It is not entirely straightforward to compare the num- 
ber density of low-ionization absorption systems at z ~ 
6 to lower redshifts. At z < 5, low-ionization sys- 
tems are generally identified by their strong H I ab- 
sorption (e.g., DLAs). A systematic search for low- 
ionization lines independent of H I column density, as 
we are forced to do here, has not been performed at 
lower redshifts. We can estimate the number den- 
sity of lower-redshift low-ionization systems, however, 
by identifying the H I column density where low- 
ionization lines, particularly O I, start to appear. O I 
is ubiquitous am ong DLAs (logiV(HI) > 20.3; e.g., 
iWolfe et al.ll2005h . and is also present in sub-DLAs down 
to logiV(H I) ~ 1 9.0 (|Dessauges-Zavadskv et~aTl 120031 : 
IPeroux et al.l 12007). Low-ionization lines are also seen 
in lower column density systems. However, these tend 
to show very strong Si IV and C IV absorption rela- 
tive to Si II and C II (e.g., iProchaska fc Buries! U999t 
iProchter et ail l2010f) . which is not seen in our systems 
(except, perhaps, for the z = 5.8865 system towards 
SDSS J1630+4012, which may contain strong Si IV). 
We therefore conclude that our z ~ 6 low-ionization 
systems are most naturally compared to lower-redshift 
DLAs and sub-DLAs with logiV(HI) > 19.0. A more 
objective comparison must await either a blind survey for 
low-ionization metal lines at lower redshifts, or a more 
complete characterization of the metal line properties of 
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Figure 12. A schematic representation of our survey results. The horizontal lines show the redshift interval surveyed for each QSO. Solid 
lines are for HIRES or MIKE data. Dotted lines are for ESI data. Filled circles represent detected low-ionization systems. The inset 
shows the cumulative probability density function (CPDF) for the survey pathlength, AX(z) (dashed line), and for the redshifts of the 
low-ionization systems (histogram). 



sub-DLAs and lower column density systems. In the fol- 
lowing sections we generally compare our z ~ 6 systems 
to lower-redshift sub-DLAs, since they are more common 
than DLAs and therefore more likely to make up the ma- 
jority of our high-redshift sample. 

The number density of DLAs has been well es- 
tablished using the large pathlen gth available from 
the Sloan Digital Sky Sur vey (|Prochaska fe Wolfel 
I2009L iNoterdaeme et alll2009f ). Surveys using higher- 
resolution spectra have meanwhile been used to deter- 
mine the number density of sub-DLAs (|Peroux et al.l 
[2005t lOTVieara et al.1120071; iGuimaraes et al.ll2009D . The 
net result is that the total number of systems with 
log A(H I) > 19.0 is 1{X) ~ 0.3 over 3 < z < 5. This is 
very similar to the number density of low-ionization sys- 
tems we find over 5.3 < z < 6.4, which suggests that 
the number density of low-ionization systems may be 
roughly constant over 3 < z < 6. This contrasts to 
the total number density of optically thick Lyman lim- 
its systems (logA(HI) > 17.2), whose number density 
per unit X increases by a factor of ^2 from z ~ 3 to 6 
(jSongaila fc Cowigl201Q[ ). It also contrasts to the num- 
ber density of highly ionized metal absorption systems 
seen in C IV, which declines by at least a factor of four 
from z ~ 3 to 6 (Paper I). 

4.1.1. Possible Redshift Evolution 

A schematic representation of the redshifts of the de- 
tected systems is shown in Figure[12j We also plot the cu- 
mulative probability density function of the survey path- 



length, AX(z) (not adjusted for completeness), and for 
the redshifts of the low-ionization systems. One peculiar 
feature is that all ten of our systems lie at z > 5.75, even 
though nearly 40% of our pathlength is at lower redshifts. 
We performed a K-S test to determine whether the dis- 
tribution of absorption redshifts is consistent with our 
pathlength-weighted survey redshift distribution func- 
tion. The test returns a D statistic of 0.39, which would 
be larger in a random sample of ten redshifts drawn from 
our survey path only 7% of the time. The distribution 
appears more unusual if we consider that, for a uniform 
distribution, the probability that all ten systems would 
randomly occur at z > 5.75 is only 0.7%. This ignores 
clustering, however. If we count only one of the systems 
towards both SDSS J1148+5251 and SDSS J0818+1722 
(i.e., seven 'groups' of systems, the maximum effect of 
clustering), then the probability increases to 3%. We 
note that completeness should not be a significant factor, 
as our completeness is somewhat better below z = 5.75 
than above. Thus, there is evidence that the number 
density of low-ionization systems may be increasing with 
redshift over z ~ 5 to 6. An expanded search at z < 5.7 
would help to clarify this. 

4.2. Mass Density 

The comoving mass density of various ions can be char- 
acterized as a fraction of the critical density at z = 0, 
Pcrit: as 

n ion = / N ion f(N ioa , X) dN ion , (3) 

C/Ocrit J 
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Table 5 

Mass Densities over 5.3 < z < 6.4 



Ion 


^ion 


[M/H] global 


Si II 


> 0.4 x 10~ 8 


> -3.9 


O I 


> 3.6 x 1CT 8 


> -3.8 


C II 


> 0.9 x 10~ 8 


> -4.1 









where f(N ion ,X) = d 2 J\f/dN ion dX . The integral in 
Eq. [3] is often approximated using the observed column 
densities within a sample as 

/ N ion f(N ion , X) dN ion ~ ^ A ^ on . (4) 

The validity of this approximation, however, is highly 
model-dependent. A power law f(Ni on ,X), for exam- 
ple, may hide most of the mass in rare, high column 
density systems. It is unclear what shape the distribu- 
tion of low ion column densities should have, as it will 
depend on the underlying hydrogen column density dis- 
tribution, the distribution of metallicities, and ionization 
corrections that are potentially important for sub-DLAs 
or somewhat weaker systems. The distribution may be 
log-normal, as there are several multiplicative factors at 
work. The present sample is too small, however, to char- 
acterize the distribution with any confidence. Given that 
much of the metal mass in the low-ionization phase may 
reside in rare systems that have both high H I column 
densities and higher than average metallicities, a con- 
servative approach is to use the present sample to set 
lower limits on f2; on using the approximation in Eq 01 
The values for Si II, O I, and C II are given in Ta- 
ble [5J We also convert these values into minimum mass- 
averaged global metallicities by dividing the ionic den- 
sities by the total mass densit y of hydrogen at z = 6, 
assuming lAsplund et al.l (|2005l ) solar abundances. No- 
tably, fHC II) is twice as high as J7(C IV) measured by 
iRvan- Weber et alJ (|2009t l. although both were derived 
from small samples of lines. The true values of f2i on may 
be factors of several higher in both cases. 

4.3. Relative Abundances 

At z ~ 6 the Lya forest is sufficiently absorbed that 
H I column densities cannot be measured for individual 
absorbers, so we are unable to directly measure metal- 
licities for these systems. We can, however, determine 
relative abundances using only the metal lines. The 
relative abundances for the systems in our high reso- 
lution sample are given in Table [51 T hese are mea- 
sured with respect to the solar values from lAsplund et all 
(200jj|), and assume no dust depletion and no ion- 
ization corrections such that, for example, [Si/O] = 
log[JV(Si ll)/N(0 I)]-log[n(Si)/n(O)j . In order to min- 
imize any potential ionization corrections, we use column 
densities measured only over the velocity interval that 
contains O I. In the case of the z — 5.8415 system to- 
wards SDSS J1623+3112, the O I and Si II A1304 lines 
appear to be blended with a C IV system at z = 4.754. 
The inverse variance-weighted mean values for the re- 
maining six systems are [Si/O] = —0.03 ± 0.03 and 
[C/O] = -0.18 ± 0.03 (la). We note that the scat- 
ter in the individual values is entirely consistent with 
the measurement errors, which suggests that the ioniza- 



tion corrections may indeed be small. T he mean val- 
ues ar e similar to the values reported in iBecker et al.1 
(2006), as expected since many of the systems in Ta- 
ble [6] were included in the earlier sample. As noted 
there, the relative abundances are consistent with enrich- 
ment from Type II supernovae of low-metallicity massiv e 
stars (IWooslev fe Weaverlll995HChieffi fc Limongil[200l . 
Our [C/O] va l ues ar e also consistent with those found by 
iPettini et aLl (120081) for metal -poor DLAs at z ~ 2 - 3 
(but see lPenprase et al.l (|2010( ). who find higher values), 
and with the [ C/O] values of low-m etallicity stars in the 
Galactic halo (jAkerman et al.l[2004l ). We will further ad- 
dress the implications of the relative abundances for stel- 
lar populations at high redshift in a forthcoming paper. 

4.4. High-Ionization Lines 

A notable feature of our low-ionization systems is the 
lack of strong absorption by high-ionization species (Si IV 
and C IV). At z ~ 6, Si IV and C IV are shifted into the 
far red and near infrared, and we stress that our spectra 
at these wavelengths are generally of modest quality. In 
most cases, however, the high-ionization lines must be 
significantly weaker than the low-ionization lines. 

Significant C IV absorption is seen for nearly all lower- 
redshift DLAs and sub-DLAs. It i s kinematically dis- 
tinct from the low-ionization lines (Wol fe fc Prochaskal 
2000), and has been interp reted as evidence for galac- 
tic winds (|Fox et al.1 [2007). We therefore investigate 
whether the non-detections of Si IV and C IV infer that 
the host galaxies of our z ~ 6 systems lack strong winds 
or other nearby, highly-ionized gas. The column densi- 
ties of high-ionization lines are plotted verses the column 
densities of low-ionization lines in Figure[13]for our z ~ 6 
systems and a sample of DL As and sub-DLAs at z < 
4.2 d r awn from the literature (iDessauges- Zavadsk v et al.1 
120031 iPeroux et al.l [20071 : iFox et al.l 120071 ). For the lit- 
erature data, direct measurements of 2V(Si II) are used 
when available. Otherwise, iV"(Si II) was estimated from 
the H I column density and metallicity measurements, 
assuming solar abundances and that all of the silicon in 
the low-ionization phase is in Si II (i.e., no ionization cor- 
rections). For the 19 sub-DLAs with Si II measurements, 
which had 19.1 < log A(H I) < 20.2, this procedure re- 
produced the published A(Si II) to within 0.03 dex in all 
cases. This is not surprising, since the metallicity esti- 
mates were often based on Si II. Column densities were 
estimated for C II using the same procedure. For one 
sub-DLA, a limit on C IV that was found to be unusu- 
ally low was re-estimated from the published spectrum. 
Otherwise, the Si IV and C IV measurements are as orig- 
inally reported. 

The systems at z ~ 6 tend to show lower high- 
ionization column densities than DLAs and sub-DLAs 
at lower redshifts. However, the upper limits on Si IV 
and C IV are consistent with a trend of declining high- 
ionization line strength with declining low-ionization col- 
umn density. This is reminiscent of the correlation 
between C IV strength and metallicity identified by 
IFox et all (|2007t ) for DLAs and sub-DLAs at z < 3.6. 
The trend suggests that the z ~ 6 systems may lack 
strong high-ionization lines because they have low metal- 
licities. Highly-ionized envelopes or winds may still be 
present but not detected in absorption with present- 
quality data. 
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Table 6 

Relative Abundances for High-Resolution Systems 



QSO 




•2-abs 


(flo,"hi) a 


[Si/O] 


[C/O] 








(km s _1 ) 






SDSS J0818- 


-1722 


5.7911 


(-86,45) 


-0.03 ± 0.05 


-0.14 ±0.04 


SDSS J0818-, 


-1722 


5.8765 


(-16,12) 


-0.11 ±0.07 


-0.16 ±0.09 


SDSS J1623- 


-3112 


5.8415 


(0,104) 


0.27 ± 1.00 


-0.30 ± 1.41 


SDSS J1148n 


-5251 


6.0115 


(-36,32) 


0.01 ±0.04 


-0.24 ±0.06 


SDSS J1148H 


-5251 


6.1312 


(-12,15) 


-0.35 ±0.35 


-0.64 ±0.40 


SDSS J1148n 


-5251 


6.1988 


(-12,10) 


-0.23 ± 0.14 


-0.33 ±0.17 


SDSS J1148- 


-5251 


6.2575 


(-10,11) 


-0.15 ± 0.15 


-0.14 ±0.20 


Mean b 








-0.03 ± 0.03 


-0.18 ±0.03 



a Minimum and maximum velocities with respect to the absorber redshift over 
which column densities were integrated 

b Inverse variance-weighted mean values excluding the z = 5.8415 system, for 
which O I and Si II A1304 appear to be blended with a lower-redshift C IV 
absorber 
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Figure 13. High-ionization (Si IV or C IV) column densities as a function of low-ionization (Si II or C II) column densities. Filled circles 
with arrows are upper limits for the z ~ 6 systems in our high resolution sample for which Si IV or C IV coverage is available. Other symbols 
are for z ~ 2 — 4 sub-DLAs (squares and large arrows) and DLAs(diamonds and small arrows) in the literature. The Si II and C II column 
densities for the sub-DLAs and DLAs are either published measurements or were estimated from published metallicities and 7V(H I) values. 
See text for details. A linear relation is plotted as a dotted line in each panel. For silicon, the relation is logiV(Si IV) = log7V(Si II) — 0.5. 
For carbon, the relation is log7V(C IV) = logA r (C II) — 1.0. The relations are not fits, but are meant to demonstrate the direction of equal 
fractional changes in the metal column densities of the low- and high-ionization phases. In this sense, the z ~ 6 systems are consistent with 
being lower-metallicity analogues of the z < 4 sub-DLAs and DLAs. We note that sub-DLAs are likely to be more common than DLAs at 
z ~ 6, just as they are at lower redshifts, and so comparing our systems with lower-redshift sub-DLAs may be the most appropriate choice. 



IRvan- Weber et al.l (pOOl and ID'Odorico et"aH (poTTf ) 
noted a possible C IV doublet at z ~ 5.79 towards 
SDSS J0818+1722, which would coincide with one of our 
low-ionization systems. There is possible C IV absorp- 
tion in our NIRSPEC data near Av = -100 km s _1 , 
where there is also Si II and C II absorption but no O I. 
The lack of O I suggests that the gas at this velocity 
may be ionized, and so some C IV might be expected. 
Confirmation must await higher quality data. 

4.5. Velocity Width Distribution 

The absorbers in our sample show a range of ve- 
locity profiles, from single, narrow lines to broad, 
multi-component systems. The distribution of Avqq 
values among the seven systems in our high res- 
olution sample is plotted in Figure Q31 We 
also show the distri bution of Atign for 28 sub- 
DLAs in the literature (Dcssaugcs-Zavad skv et al. l l2003l ; 



IPeroux et aD l2007t iFox et al l I20071) . For the systems 
in iDessauges-Zavadskv et al.l (2003) and IPeroux et al.1 



(2007), the velocity widths were calculated from pub- 
lished Voigt profile fits. There is some indication that 
the low-ionization systems at z ~ 6 are more likely to 
be narrow. A K-S test, however, indicates that there is 
a reasonable likelihood (27%) that the two samples were 
drawn from the same parent distribution. 

4.6. Mass-Metallicity Relation 

The high opacity of the Lya forest at z ~ 6 prevents 
us from measuring H I column densities for these sys- 
tems. We are therefore unable to obtain metallicities di- 
rectly. Previous works have noted correla tions between 
metallicity and either the velocity width ([Ledoux et al.l 
20061) or the rest equiv alent width of Si II A1526 (W 152 6] 
Prochaska et al.l 120081) . These relationships have only 
been calibrated for DLAs, however, while our systems 
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Figure 14. Distribution of Avgo values for the seven z ~ 6 
systems with high-resolution data (shaded histogram) and for 
lower-redshift sub-DLAs (filled histogra m). The 28 sub-DLA mea- 
surem ents were computed from data in Dcssaugcs-Zavadsky ct al. 
(2003) and IPeroux et al.1 l[2007t) . There is some indication that 
narrow systems are more common at z ~ 6. Both samples are 
small, however, and the overall distributions are statistically not 
significantly different. 

potentially have lower H I column densities. Some in- 
sight can be gained, however, by comparing the relation- 
ship between Aw go and Wi^e m our study to systems at 
lower redshifts. 

Figure [15] shows the relationship between Ai>g and 
W1526 for DLAs and sub-DLAs at z < 4.5 and for 
our low- i onization system s at z ~ 6. The DLA data 
are from iProchaska et al.l (|2008|) . Values of Avg and 
W1526 for sub-DLAs were computed from pu b lished 
Voig t profile fits f r om iDessauges -Zavads kv et al.l (2003) 
and IPeroux et al.l (|2007l ). We also include the three 
metal-poor DLAs a nd one metal-poor sub-DLA from 
iPettini et al.l (|2008| ). The z ~ 6 sample includes the 
seven syste ms in our high res ol ution sample. 

Following IProchaska et al.l ()2008l) . we plot the ex- 
pected relation for fully saturated absorption, VFff^ = 
(1526.7 A) Av 90 /c, as a dashed line in each panel. We 
also show this relation scaled by a factor such that mean 
residual in log PF1526 is zero for each data set. The DLAs 
roughly follow the relation expected for saturated absorp- 
tion, while VF1526 for the sub-DLAs is lower by 0.44 dex. 
For the z ~ 6 systems, the offset is 0.81 dex, or a fac- 
tor of ~2 lower than the sub-DLAs. At a given velocity 
width, therefore, the low-ionization systems have signifi- 
cantly weaker metal lines than either DLAs or sub-DLAs 
at z < 4.5. 

We can translate the offset in WV^fi into an appro xi- 
mate difference in metallicity. IProchaska et al.1 (|2008l ) fit 
a relation for DLAs of the form [M/H] = a + b log VF1526, 
finding (a, 6)dla = (—0.92,1.41) with an r.m.s. scat- 
ter of 0.25 dex. We have fit a similar relation for 
the sub-DLAs using the published metallicities, finding 
(a, 6) sub -DLA = (-0.28 ± 0.15, 1.19 ± 0.18), with a scatter 
of 0.33 dex. It is apparent, therefore, that the metallicity 
estimate will depend on the H I column density of the 
absorber. If we assume that our z ~ 6 systems are at 
least sub-DLAs, then the decline in W1526 suggests that 
the metallicity of the z ~ 6 is systems at a given Ai>go is 
lower by >0.4 dex than that of lower-redshift systems. 



The velocity width is expected to correlate with the 
mass of the dark matter halo (e.g. . lHaehnelt et al.lll998t 
Mailer et ail 120011: IProchaska et al.l 120081: iPontzen et al l 
2008). This decline in W1526 at a given Avgo can there- 
fore be interpreted as evolution in the mass-metallicity 
relation for the host galaxies of low-ionization absorbers. 
iLedoux et al.l (|2006l ) noted a similar redshift evolution 
for DLAs in their [M/Hj-Awgo relation, in which the 
metallicity at a given Awgo declined by ~0.3 dex from 
1.7 < z < 2.43 to 2.43 < z < 4.3. Our results demon- 
strate that this decline in the mass-metallicity relation 
likely extends to galaxies near the reionization epoch. 

5. DISCUSSION 

We now consider whether the number density of 
low-ionization systems is expected to remain roughly 
constant with redshift, as observed. Previous works 
have shown that the number density (per unit X) 
of both DLAs an d sub-DLAs shows little evolution 
over 3 < z < 5 (IPeroux et all 120051: IProchaska et al] 
20051: IProchaska fc Wolfdl2009t iNoterdaeme et al.l 120091 : 
Guimaraes et al. 2009). Our observations suggests that 
the number density remains roughly constant out to 
z ~ 6, but why should this happen? The number den- 
sity of a population of absorbers will depend on their 
comoving number density, <p c , and their mean physical 
absorption cross section, (cr a bs), as 

£{X) = ^-4>c (a ab s) • (5) 
tip 

Observation al (e.g.. iCooke et all l2006alfbl) and theo- 
retical (e.g.. iGardner e t al.l Il997t l Haehnelt et ail 119981: 
iNagamine et alll2004D arguments generally suggest that 
DLAs are hosted by dark matter haloes of mass M > 
10 8-9 M Q , and may be largely as sociated with halos i n 
the range 10 9 < M/M Q < 10 11 (jPontzen et all 120081 
For a Sheth-Tormen halo mass function, 4> c for 10 10 M 
halos, for example, will be lower at z = 6 than at z = 3 
by a factor of ~4 (ISheth et all 120011: IReed et al.l 12001 . 
It is possible that the host halos of low-ionization sys- 
tems are somewhat less massive at higher redshifts, which 
would increase their space density. The similarity of the 
velocity width distributions of the z ~ 6 systems and 
low-redshift sub-DLAs suggests that the halo masses are 
similar, although the velocity width may only be mod- 
erately se nsitive to the mass. For an NFW dark mat- 
ter profile ([Navarro et al.lll~997l ). the virial velocity scales 
with the halo mass as iVir oc M v / 3 iJ(z) 1 / 3 , where H(z) is 
the Hubble parameter at redshift z. If At;gn scales with 
the virial velocity (e.g.. lHaehnelt et aLlll998j ). then a de- 
crease in the typical halo mass by a factor of ten from 
z = 3 to 6 would only produce a decline in the typical 
Augo of ~0.2 dex. Such a decline is consistent with the 
present data, but as noted in Section l4.5l the distribution 
of Augo values for the z ~ 6 sample is formally consistent 
with that of lower-redshift sub-DLAs. A larger sample is 
needed to determine whether there is a genuine decrease 
in the mean velocity width. 

Multiple factors may increase the absorption cross- 
section with redshift. At a given mass, an NFW profile 
will have a higher density, particularly in the inner re- 
gions. The e xternal UV background also decreases from 
z ~ 3 to 6 dFan et al l 1 20061: iBolton fc Haehneitl l2007t 
iWvithe fc Boltonl 120101: ICalverlev et all 120101 ) . Both of 
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Figure 15. Rest equivalent width of Si II A1526 as a function of Ango. The left an d middle panels are for z < 4.5 DLAs and z < 4.2 
sub-DLAs, respectively, from the literature. These include the metal-poor systems from Pcttini et al. (2008), which are shown as diamonds. 
See text for details. The right panel is for our z ~ 6 systems with high-resolution data. The dashed line in each panel shows the relation 
expected for fully saturated absorption, Wfjjjg = (1526.7 A)Ati9o/c. The dotted lines show this relation offset vertically such that the mean 
deviation in log W1526 is zero. For the DLAs the offset is close to zero. At a given velocity width, the metal lines for the z ~ 6 systems are 
significantly weaker than for either the lower-redshift DLAs or sub-DLAs. 



these factors should increase the radius out to which a 
halo of a given mass can become self-shielded. 

We explored simple models using NFW profiles to pre- 
dict how the cross-section of low-ionization gas should 
evolve with redshift. Assuming that the gas traces the 
dark matter and is in photoionization equilibrium, it is 
straightforward to predict the radius at which the gas be- 
comes self-shielded and to estimate the cross-section over 
which the projected H I column density exceeds a given 
value. For this we assumed a gas temperature of 20,000 K 
and a metagalactic H I ionization rate T = 10~ 119 s _1 
(10~ 12 - 8 s" 1 ) at z = 3 (6). This simple model suggests 
that the cross section over which a halo of a given mass 
would be a DLA should scale roughly linearly with the 
halo mass, which agrees broadly with recent theoreti- 
cal models (iNagamine et al.l 120041: iPontzen et all 120081: 
iTescari et al.l 120091 ) . It also predicts that the DLA cross 
section for a given mass halo should increase by a factor 
of ~4 from z = 3 to 6. If DLAs are mainly hosted by 
dark matter halos with masses M > 10 9 M , then this 
increase in the cross section will offset the decline in the 
comoving spatial density to produce a roughly constant 
number density of DLAs with redshift. The same is true 
for sub-DLAs in this model. We stress that this is an 
extremely simple model, and it is unlikely to reproduce 
detailed observations such as the column density or ve- 
locity width distributions of DLAs. Nevertheless, it gives 
some assurance that the apparent lack of evolution in the 
number density of low-ionization absorbers can be plau- 
sibly explained by a decline in the spatial density of dark 
matter halos combined with an increase in the mean ab- 
sorption cross section. More detailed analysis must be 
left for future work. 

The high number density of low-ionization systems 
makes it likely that they arise from galaxies below the 
detection limits of current galaxy surveys. Integrating 
the z ~ 6 luminosity f unction of i-dropout galaxies from 
iBouwens et al.1 (|2007l ) down to the limit of the Hub- 
ble Ultra Deep Field (~0.04L* =3 ) gives a space density 
4>(Mab < — 18) = 5 x 10~ 3 per comoving Mpc 3 . For 
£(X) ~ 0.25, this implies a mean absorption cross-section 



(f7 a bs) ~ 1 x 10 4 kpc~ 2 , or (i?abs) — 60 kpc for a circular 
projected geometry. The Lya-emitting galaxies (LAEs) 
at z ~ 6 included in current surv e ys wo uld need an even 
larger cross-section. lOuchi et al.l (|2008f) found a density 
of LAEs at z = 5.7 of ~7 x 10~ 4 per Mpc -3 down to a 
Lya luminosity of 2.5 x 10 42 erg s _1 . These LAEs would 
require (i? a bs) — 160 kpc to account for the observed 
number density of low-ionization absorption systems. In 
contrast, constraints from paired QSO lines of sight sug- 
;est a typical DLA size at z = 1-2 of -Rdla — 5 — 10 kpc 
Monier et aLl200alCooke et alj2010D, although in some 
cases they may be larger (e.g.. iBriggs et al.lll989f ). The 
sizes of sub-DLAs are less well constrained, but they may 
plausibly be a factor of two lar ger than DLAs (see , for 
example, th e scaling relation in|Monier et all (|2009f ). or 
Figure 2 in IPontzen et ail (|2008l) h In that case, they 
would still need to be significantly larger at z = 6 to 
arise solely from galaxies similar to those in current sur- 
veys. It seems more plausible, therefore, that these low- 
ionization systems trace fainter, more nume rous sources. 
This w ould be consistent with the results of lRauch et al.1 
(2008), who find a likely connection between DLAs and 
a population of ultra- faint LAEs at z ~ 3. 

We note that these absorption systems are poten- 
tially our first direct probe of the 'typical' galaxies that 
are responsible for hydrogen reionization. It is widely 
believed that in order for reionization to complete by 
z = 6, the majority of ionizing photons must come 
from galaxies below the limit s of current imaging sur- 
veys (e.g., IBouwens et al.ll2007t iBolton fc Haehneltil 20071 
20091 iBunker et al l 120101: iMcLure et al l 
120101 ) . Metal absorption lines may 
therefore be the most efficient means of detecting these 
sources, as directly detecting them in emission must 
await facilities such as ALMA and JWST. This study 
already delivers insights into the properties of reion- 
ization galaxies, including the suggestion that they are 
likely to be more metal poor than their lower-redshift 
counterparts, at least in the gas phase. If their stars 
also have low metallicities, then they may be highly 
efficient in producing ionizing photons (jSchaererl 120031 : 
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Venkatesan fc Truranl [2003). As noted by others (e.g., 
Ouchi et al.ll2009ft . this would potentially help to explain 
how reionization is completed by z ~ 6. We add that, 
in addition to their role in reionization, these low-mass 
galaxies may contribute si gnificantly to the c hemical en- 
richment of the IGM (e.g.. Booth et al.ll2010h . 

6. SUMMARY 

We have conducted a search for low-ionization metal 
absorption systems spanning 5.3 < z < 6.4. Our sur- 
vey includes high- and moderate-resolution spectra of 
17 QSOs with emission redshifts z cm — 5.8 — 6.4. The 
total survey pathlength is AX — 39.5 (Az = 8.0, or 
Al = 3.6 comoving Gpc), of which roughly half is covered 
by high-resolution data. We searched for low-ionization 
systems by looking for coincidences in redshift between 
Si II, C II, and O I. In total we dete ct ten systems , five o f 
which were previously reported bv iBecker et al.l ((2006). 
Each contains Si II and C II, and all but one contains 
O I. The majority are detected in our high-resolution 
data, which is consistent with the fact that these data 
are more sensitive to narrow absorption lines. 

The line-of-sight number density of absorption sys- 
tems, uncorrected for completeness, is £(X) = 0.25t ^i3 
(95% confidence). This is similar to the number density 
over 3 < z < 5 of DLAs and sub-DLAs (log N(R I) > 
19.0), which constitute the main population of low- 
ionization absorbers at those redshifts. The fact that 
the number density of low-ionization absorbers is roughly 
constant out to z ~ 6 is in sharp contrast with the 
evolution of high-ionization systems traced by C IV, 
which show a marked decline at z > 5.3 (Paper I; 
iRvan- Weber et all 120091) . At z ~ 6, low-ionization sys- 
tems with log N(C II) > 13 are more common than high- 
ionization systems with logiV(ClV) > 13, a reversal 
from lower redshifts. 

The roughly constant number density of low-ionization 
systems over 3 < z < 6 may be explained if they are 
hosted by lower-mass dark matter halos at higher red- 
shifts. Alternatively, if the systems at z ~ 6 are hosted 
by halos with masses similar to those that host DLAs 
and sub-DLAs at lower redshifts, then the apparent lack 
of evolution may occur if dark matter halos of a given 
mass have a larger mean cross-section of low-ionization 
gas at higher redshifts due to the higher gas densities and 
weaker UV background. 

Although the H I column densities cannot be measured 
for these systems, we are able to infer some of their prop- 
erties by comparing the velocity widths and metal line 
strengths to samples at lower redshifts. For the seven 
systems with high resolution data, the velocity widths 
span a similar range as sub-DLAs at 2 < z < 4, although 
there is some indication that the z ~ 6 systems tend to 
be narrower. The lines in the z ~ 6 systems are also 
weaker in the sense that, at a given velocity width, the 
inferred equivalent width of Si II A1526 is lower than 
for z < 4 sub-DLAs by a factor of two, and by a fac- 
tor of six compared to z < 4 DLAs. This implies that 
the mass-metallicity relation of the host galaxies evolves 
towards lower metallicities at higher reds hifts, a trend 
that has also been noted over 2 < z < 4 (|Ledoux et al.l 
2006). Assuming the z ~ 6 systems span a similar range 
in log N(H I) as the low-ionization systems at z < 4, the 
strength of the metal lines implies a decline in the gas- 



6 of at least 



phase metallicity from 2 < z < 4 to z 
-0.4 dex. 

The mean relative abundances are [Si/O] = — 0.03 ± 
0.03 and [C/O] = -0.18 ± 0.03 (ler uncertainties), 
assuming no dust depletion and no ionization correc- 
tions. These are co nsistent with t he va lues found for 
metal-poor DLAs by iPettini et all (120081) and fo r low- 
metallicit y halo stars by Akerman et all ([2004). As 
noted bv Becker etal] (|2006h . the relative abundances 
are broadly consistent with enrichment from Type II su- 
pernovae of low-metallicity massive stars. 

Our z — 6 systems are also notable in that they lack 
strong high-ionization lines (Si IV and C IV), which are 
ubiquitous amon g lower-redshift DLAs and sub-DLAs 
(|Fox et al.l 120071) . The absence of these lines is consis- 
tent with a similar fractional decline in the metallic- 
ity of the low- and high-ionization phases, however, and 
does not necessarily indicate that these systems lack ha- 
los of highly-ionized gas. Deeper spectra will help to 
determine whether the z — 6 absorbers have signifi- 
cantly lower N(Si IV)/iV(Si II) and N(C IV)/A(C II) ra- 
tios than lower-redshift systems. 

The overall consistency between the properties of the 
z — 6 systems and those of lower-redshift DLAs and sub- 
DLAs suggests that the z — 6 absorbers arise from galaxy 
halos, rather than from remnants of neutral IGM at the 
tail end of hydrogen reionization. One notable aspect of 
our survey, however, is that all ten of our systems occur 
at z > 5.75, while roughly 40% of our pathlength is at 
lower redshifts. Expanded surveys over 4.5 < z < 5.7 
will help to determine whether the number density of 
low-ionization systems truly remains constant out to z — 
6, or whether there is a decline up to z — 5.7 followed 
by an increase at higher redshifts. The latter scenario 
would potentially indicate a strong evolution in the UV 
background near z — 6. 

Finally, the high-number density of low-ionization sys- 
tems at z ~ 6 suggests that we are probing galaxies be- 
low the detection limits of current i-dropout and Lya- 
emission galaxy surveys. As such, these absorption sys- 
tems are potentially the first observations of 'typical' 
galaxies responsible for hydrogen reionization. The low 
metallicities we infer suggest that these galaxies may be 
highly efficient at producing ionizing radiation, a fact 
which would help to explain how the IGM becomes fully 
ionized by z — 6. As more and higher- redshift QSOs 
are discovered, absorption lines will continue to provide a 
unique probe of the reionization era that will complement 
studies with ALMA, JWST, and other next-generation 
facilities. 
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